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Nomenclature 

» burning rate . 

= pressure-, temperature-sensitivity coefficient 
= coefficients in perturbation solution, where i 
is equal to 1, 2, .... 10 
= rate -of- strain tensor 

= Froude number ... . 

B inverse Froude number (gravitational acceleration) 
= perturbation wave number 
= unit normal 
= Prandtl number 
= pressure 

= quantity defined following Eq. (26) 

= quantity defined following Eq. (26) 

35 time variable 
s= perturbation velocity vector 
s velocity vector 
= moving coordinate system 
= surface-tension coefficient 
s small bookkeeping parameter . . i 

s perturbation pressure 
as gas-to-liquid thermal diffusivity ratio 
= gas-to-liquid viscosity ratio 
= gas-to-liquid density ratio 

= location of gas/liquid interface 

b perturbation in location of gas/liquid interface 
= complex perturbation frequency 


Subscripts and Superscripts 

i = inner wave number regime or integer variable 


wl 1/ = far outer wav^umber regime 

l 
g 


/ = liquid 

g = gas 

i l._aV o = outer wav£)umber regime 

*1^" * = scaled quantity 
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Introduction 

YDRODYNAM1C (Landau 1 ) instability in combustion is I typ- 
jl jl ically associated with the onset of wrinkling of a flame surface, 
which corresponds lo the formation of steady cellule structures as 
the stability threshold is crossed. This type of instability was ong- 
inaily described by Landau 1 and is attributed to thermal expansion 
across a combustion front. Although gaseous combustion was de- 
termined lo be intrinsically unstable in Landau s analysis, u was 
demonstrated that additional effects, such as gravity and surface 
tension, that enter when the unbumed mixture is a liquid result in n 
specific stability criterion. However, dus analysis, along with a later 
study by Lcvich 2 that considered viscous effects in lieu .of surface 
tension, assumed that the the combustion front propagated normal to 
itself with constant speed, whereas it is now recognized that there 
a dynamic interaction between the burning rate and local conditions 

31 Fortoc problems in wliicli pyrolysis, cxodicrmic decomposi- 
tion and/or combustion occurs in an intrusive region in the vicinity 
of the liquid/gas interface, the dynamical coupling of die burning 
rate with the underlying hydrodynamics of the flow can be achieved 
through an analysis of the thin combustion/interface region. An al- 
ternative approach, however, is to simply postulate a phenomeno- 
logical dependence of the local burning rate on pressure and tem- 
perature and to obtain results in terms of suitably defined sensitiv- 
ity parameters. Both types of mediodologies have been . applied Uo 
die problem of solid-propellant combustion, and each offers certain 
advantages. 1,4 In die present scries of studies on hquid-propcllaiU 
combustion, 1 -’ die latter approach has been adopted, thereby gen- 
eralizing the Landau'/Levich 2 model to allow for a coupling of die 
burning rate with the localpressure and temperature liefc ^ . 

af the Summarizing Jrfsome of the results obtained from the . 
present model, it has been shown dial when only the pressure sen- 
sitivity of the burning rale is taken into account an appropriately 
generalized stability witerion for cellular (Landau 1 ) instability is ob- 
fained Exploiting the realistic limit of small gas-to-liquid density ., 
ratios it isfeund that the stable region occurs.for negative values of 
the pressure sensitivity parameter, with the origmal L^tdau modd 
being intrinsically unstable in this lnmt. In particular, die neutral 
stability boundary possesses a local minimum when plotted against . 
die disluibancc wave number, winch suggests that as the pressure- 
sensitivity parameter decreases in magnitude, the hqu.^gas iiUcr- 
face/front develops cells corresponding to classical hydrodynamic 

This report is a preprint of an article submitted to a journal for 
publication. Because of changes that may be made before formal 
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instability. 5 This minimum rcflcctsJhat surface tension and viscos- 
ity are stabilizing influences for slftrt-wave disturbances, whereas 
gravity is a stabilizing influence for long-wave perturbations. As a 
result, the effect of reducing the gravitational acceleration to rru- 
crogravity levels is to shift the neutral stability minimum to smaUer 
wave numbers. Thus, in the microgravity regime, Landau instability 
becomes a long-wave instability phenomenon, implying the appear- 
ance of large cells along the combustion interface. ; 

Aside from the classical cellular form of hydrodynamic instabil- 
ity this dynamic generalization of the Landau /Levich model also 
predicts the appearance of a pulsating form of hydrodynamic in- 
stability, corresponding to the onset of temporal oscillations in the 
location of the liquid/gas interface. This form of hydrodynamic in- 
stability occurs for negative values of the pressure-sensitivity param- 
eter that are sufficiently large in magnitude. 6 Consequently, stable, 
planar burning is predicted to occur in a range of negative pressure 
sensitivities that lies below the cellular boundary and above the pul- 
sating boundary just described. A stable range of negative pressure 
sensitivities is applicable, for example, to certain types of hydroxy- 
lammonium nitrate- (HAN-) based liquid propeflantsforwhich non- 
steady modes of combustion have been observed. 8 -* (Though less 
common, ranges of negative overall reaction orders/pressure sensi- 
tivities have been reported for sufflciently diluted gaseous hydrocar- 
bon fuels as well 10 ). The appearance of both pulsating and cellular 
forms of hydrodynamic instability is analogous to the two corre- 
sponding types of thermal/diffusive instabilities that occur for suf- 
ficiently large and sufficiently small Lewis numbers, respectively. 

When the effect of a temperature sensitivity in the burning rate 
is included in the analysis, substantial modifications to the preced- 
ing stability description can occur. Specifically, if die temperature- 
sensitivity parameter is sufficiently large relative to the parameter 
corresponding to pressure sensitivity, the pulsating hydrodynamic 
stability boundary can develop a turning point, that is, become 
C shaped, in the (disturbance- wav£* umber, pressure-sensitivity) 
plane. In that case, the stable region for small wave numbers dis- 
appears, and liquid- propellant combustion is predicted to be intrin- 
sically unstable to the nonsteady form of hydrodynamic instability 
for all sufficiently large disturbance wavelengths. This has been de- 
scribed in detail in the limit of zero viscosity, 7 and the purpose of the 
present work is to extend that analysis to the fully viscous model. 
Viscous effects were shown to have a substantial influence m the 
absence of thermal sensitivity, where it turned out that the stable 
region became significantly widened when viscosity was present, 
and the same result will be demonstrated when diermal effects are 
present. However, the same intrinsic pulsating instability that occurs 
for sufficiently large temperature sensitivities and sufficiently small 
wave numbers in the inviscid case will be shown to be preserved 
even when viscosity is included. These results lend further support 
to the notion that a likely form of hydrodynamic instability in hquid- 
propellant combustion is of a nonsteady, long-wave nature, distinct 
from the steady, cellular form originally predicted by Landau. 

The physical nature of the pulsating form of hydrodynamic m- 
stabihty described here, like the pulsating form of therm al/diffiisive 
(or reactive/diffusive) instability: is manifested through an oscil- 
latory imbalance between reaction-front perturbations and those 
processes that act to dampen such perturbations. In die case of a 
pulsating reactive/diffusive instability, such as occurs in gaseous 
combustion, 11 smaller mass-lo- thermal diffusivity ratios, that is, 

„ larger Lewis numbers, allow a relatively greater concentration ot 
reactant in the reaction zone. This in turn triggers, for sufficient y 
; large Zel'dovich numbers, a more intense reacuon, which leads to 
an imbalance between temperature perturbations that accelerate the 
front and cause the profiles to steepen, and diffusion, which trans- 
V f crs heat to the unbumed mixture and thereby reduces the reaction 
■■ 1 intensity. In the purely hydrodynamic problem, a negative pressure 
' sensitivity plays a somewhat analogous role to that of diffusion be- 
cause positive pressure perturbations will either locally accelerate 
or decelerate the front, depending on whether the pressure sensi- 
v tivity is positive or negative. Thus, positive pressure sensitivities 
lead to intrinsic instability, 5 whereas negative pressure sensitivities 
that are sufficiently large in magnitude lead to an overcorrecUon 
in the local burning rate in response to a hydrodynamic pressure 


disturbance. In the latter case, an oscillatory imbalance between hy- 
drodynamic perturbations and corresponding variations m the local 
burning rate is/thus/established. As indicated by the subsequent f 
results, flic inciusiorf of viscosity and a thermal sensitivity in the re- 
action rate, where the latter results in a coupling of the thermal and 
hydrodynamic fields, accentuates this effect through the inclusion 
of thermal/diffusive processes as already described. 

Summary of the Mathematical Model 
The mathematical model was described previously, 5,17 but is 
briefly summarized here for completeness. Specifically, it is as- 
sumed that the combustion front coincides with the liquid/gas inter- 
face, where pyrolysis and/or exothermic decomposition occurs. De- 
noting the nondimensional location of this downward-propagating 
interface by x 3 = <K*i, * 2 , 0,where* 3 is the vertical coordinate and 
the adopted coordinate system is fixed with respect to the stationary 
liouid at jci = —oo, we transform to the moving coordinate system 

x=x lt y = x 2 , and z = x 3 - such Uut lic * uid/g “ 

interface always lies at z = 0. Conservation of mass, energy, and 
momentum within each phase then gives 


V • v = 0, 


z*0 


£® _!±£® +K . V 0 = ( ‘ )v J 0 . 
Bt Bi Bi l A J 

_ii^ + (K.V)^ = (0.0. -Fr' 1 ) 
Bt dt Bz 


l 

P _l . 


z£0 
Vp 




zSO 


( 1 ) 

( 2 ) 


(3) 


where 0 is temperature, l’n and Pr t are die liquid- and gas- 
phase Prandd numbers, X and c (to be used) are thermal diffusivity 
and heat-capacity rados, and Fr' is the inverse Froude number 
(eravitational acceleration). Other nondimensional parameters in- 
troduced subsequently include the gas-to-liquid viscosity ratio p 
(n\Pr = aPr,) and the unbumed-to-bumed temperature ratio cr„. 

Equations (1-3) are subject to a set of boundary and interface 
conditions given by 


y = o, 0 = 0 at 

© = 1 at z = +oo, 
H, X V- = «, x F+ 


Z = -oo 

0U-O- =01,-0. 
3<J> 


(4) 

(5) 

( 6 ) 
( 7 ) 


n, ■ (»>_ - pv+) = (1 - P)SW-gf 

n. ■ v. - S(<t>)^r = A(p ©It-o) / , 
at 

p\i- 0 - - pl,-o* = ", • l PM", • v +> - •'-("« • *’-) 

34> / 

- pXPr,e + ■ ii, + Pne- • ii,] + • (f- ~ pv+)S(<J>) — 


-y(-Vii,) 


r 

ii, x pM«, • *+) - - *’-> + (v - - 

= h, x (pkPr,e+ • ft, — Pne- • ft,) 
h, • (cpXV 0|j«o+ “ V©| t .o-) 

= h, • l(cpv+ - v-)©I t -o + c(o u pv+ - v-)] 

3<J> 


-f [(I — cp)©| t »o + c(l — cr B p)]5(4>) 


dt ' 

. .Id. i 


e = 


( 8 ) 


(9) 


.1 - <* U f; 
( 10 ) 


’ ' d*T* 


where = v\ t -o* . = * I* - o* * ] &P- t 5 ' 1 0) c ?™' 5pond t0 c ° 

linuity of tlie Lransverse velocity components (no-shp). conservation 
of (normal) mass flux, the mass burning rate (pyrolysis) law, 
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conservation of flux of the normal and transverse components 
of momentum, and conservation of heat flux. Here, 5(<P) - 
n + 4 . <|> 2 )" 1/f2 , the unit normal h s = (— ^ — 

the expressions for the gradient operator, the Laplacian, and 
curvature in the moving coordinate system are given by V- 

(a, - a, - <Mz, a«). vl = + V- (1 + # t, V& I 

2<b,3„ - 2<b,3„ - (4>„ + 4>„)3,. and -V • n, = <*>«(! + <t>,) + 
<j> (1 + 4 > 2 ) - 2<J> J 4> y <S>x y . respectively. However, the vector v still 
denotes the* velocity with respect to the (xj , x 2 , x 3 ) coordinate sys- 
tem Finally, we note that the nondimensional mass burning rate 
appearing in Eq. (7) is assumed to be functionally dependent on 
both the local pressure and temperature at the liquid/gas interface. 
By definition, A = 1 for the case of steady, planar burning, but per- 
turbations in pressure and/or temperature result in corresponding 
perturbations in the local mass burning rate. 

Because the thermal and hydrodynamic fields are coupled only 
through the temperature dependence of the mass burning rate A ap- 
pearing in Eq. (7), the strictly hydrodynamic problem for p, v, and 
<p. can be analyzed separately when A is assumed to depend on 
pressure only. 5,1 * In the present work, wc focus on the fully coupled 
problem to determine how the hydrodynamic stability boundaries 
are modified when the local burning rate depends on temperature as 
well as pressure. Our stability results wiilyihus/depend on two sensi- 
tivity parameters, A p and A©, defined as A p = dA/dp |e « 

Ae = 3 A/ 30 | 0 -i.p-o. where 0*1. and p=0 are the interface 
values of temperature and pressure of the basic solution in Eq. (I )■ 
Though an explicit expression for the reaction rate A is not needed 
in the present analysis, we note that, because the nondimensional 
activation energy is typically large, the temperature sensitivity A e 
would likely be larger in magnitude than the pressure sensitivity 
A , which will have some bearing on the relative scalings of these 
parameters that will emerge in the following analysis. 

Basic Solution and Its Linear Stability 
The nontrivial basic solution of the preceding problem that corre- 
sponds to the special case of a steady, planar deflagration is given by 

= -t, v° = (0, 0, v°) 


z < 0 

©°(z) = 

fe z , 

o 

V 

N 

r l - 1, z>0, 

It. 

z > 0 

f -Fr'z + p-' 

- 1, 

z<0 


P {z)= \-pFr-'z, 


z > 0 

(11) 


wi| 2 »=o- “ Uilz=»o+ — iP 1 
Milt-o- -Hilt-o* = CP -1 ~ 1 )4>y 

W 3 1* = 0“ - P«3lx-0+ = (1 - P)h 
Mj| x .o- -0f = ApfU-o* + A©0U,.o+ 

£ | t ,0 CU-0+ = 2CM3U-C+ ~ Will-0*) 


/ Buy 9 m 3 

+2Pr \ir 


- 2(1 - p)<t>, 


- y(<Pxx +^jr) 


„0 + - x-0 + 


/ du\ 

“Mir 


( Bu 2 I , 3“3 

irL + i7 t . 0 . 


The problem governing its linear stability may be formulated, be- 
fore introducing any further approximations, in a standard fashion 
by introducing the perturbation quantities y,/) = < l > (x,)',z,0 
' — <S>°(/) «(x,y,z,/) = v(x,y,z,t ) — y°U). £(x,y,z,0 = 

pU,y,z./)-p°(z). and 0{x,y, z, <) = © - ©“(a ) -p /&• 
The problem obtained when Eqs. (1-10) are hneanzed about the 
basic solution (11) is then given in terms of these perturbation van- 
ables by 

3^ £«, z,40 U2) 

dx dy dz 

■ r i i a* du (h l 1 l Fr -i£i/^Tl 

Jy / ip I v>* M axy yj 
^ t /ac r^T7a* *<YJ Pri U* + 3i “ + 

V 1 p r T *> t ^ 1^ ^ 3z! b 

^ z§ o (13) 

'^^fil33 ae i !/£!£, £!*+!!£) 

{ p j a7 + 3z ” I. o ) + \p 4V 3 * 2 3 > J 3z Vftf) 

zSO (14) 

u = 0, 0 = 0 at z = -oo 

0=0 at z = +oo, B\t-o+ — 3 li»o- = 4>i (15) 


C pX — I ———I — c3|,»o+ + 3|i-o U3li-o-+^< 

■ (21) 

where Eqs. (16) and (17) have been used to simplify Eqs. (18-21). 

Nontrivial harmonic solutions for 4>. u. and f, propor- 
tional to exp(io>r +ik l x + ik 2 y), that satisfy Eqs. (12-14) and 
the boundary/boundedness conditions at z = ±°o are given by 
^ = exp(iai( +ikix+ikiy) and - , 

\b t t* z -Fr', ■' K 0 

C = exp(i<u/ + i*i.t + ‘*iy)L e -^-pFr~\ Z>0 ( 22 ) 

uj = cxpO'arf + ik^x + ik 2 y) i 

■ (bje* 1 -ikdico + ky'bie^, z<0 

x - /Jt,(i<up -ky'bie-H, z > 0 (23) 

u 2 = exp(i<o/ + ikix +ikiy) 

(b 5 e« l - ikiiioi + ky'bie 11 , i< 0 j 

x | he 11 - iki(io>p - *)'%«“*’. , z > 0 (24) 

hj = exp(i<u( + ik)X + ikty) • ! 

[he** -Hiw + ky'b^ 1 , z< 0 
X \b,e rl + Hieop - ky'bie ~ u , z>0 (25) 

; i 

$ = exp(io)/ + ik\x + ik 2 y) 

• - [iai + Jl 2 -q{q + l)]'‘^ eto . + ln , 

x +kWoi) 1 -k 1 ]-'b,e' l+ 'K : z<0 ^ 

b lo e“, z>0 ' 

, ft' ? : - .r 

where k = {k] + k])' 11 is die overall disturbance' wave 
and the quantities p. q, r. and s are defined as • P— 

4 Uto + k 1 )] ,fl , 2/+/ <7 = 1 + [1 + 4Pr,(iaj + /’r, k 2 )] ul , 2 (tPr,r ~ : 
1 - u + 4 vPntui + nPnk')]'» and ; 2 pks = l - U +; 
4 + X* J )1 1/2 . Substituting this solution into the interface con- 

ditions (16-21) and using Eq. (12) for z =* 0 y ield 1 1 commons for 
die 10 coefficients b,-b, 0 and tire complex frequency (dispersion 
relation) iio(k). In particular, these conditions .are given by 
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ik\b-$ + ikibs +qbj = i'fti&4 + ik 2 b$ + r ^8 — 0 
£?3 “ z". fri — &4 + ■ 


ico + A: 


— k 


b s - bt + t— — - fci = (P~‘ - l)i*J 

5 jo> + * top-* 

bj — T—rrbi — pb 2 — -r~ — = (1 — P) ,a> 


(27) 


(28) 


to + A: 


/ cd/? — £ 


X 1 

b, -i,, - A p bi - A© £?io =io>- pFr A p 

. ia) + k 

r, + “211,,, _ f| + Ifr, - 2Pr,qbj 

[_ ia> + A;J [_ ia)p — k J 

- 2(1 - p - tiPr,r)b, = (1 - pXFr' 1 - 2io) + yk 1 

uPorbi + ^r^—p- bi + ik\nPribt - Pnqby + ■ - — -rb { 


(29) 


(30) 


iwp - k 

— IkiPrt 67 = 0 

^rb^Jh^ 


id) + A: 


(31) 


likjkPn 

b 2 + ik 2 pPrib% - /Vtffe + + 


i<t)p - * 

-ik 2 Pr,b 1 = 0 ( 32 ) 

b,„ - 1, + [/o + * 3 - ?(? + Dry - w»? - *r% = 1 

(33) 

-7 + 1 


V-c + c P Xs)b l0 -pb,+ iaj + Jk2 _ <?((?+1) 


+ l\bi 


ico + k 




1 + to 


(34) 


// 


34). These first emerged in our analysis of cellular instability using 
the generalized model in the limit A© = 0, but they are also relevant 
when one considers the pulsating form of instability and when A© is 
allowed to be nonzero. In particular, in the case of cellular instability 
and zero temperature sensitivity, there are three wav^mmber scales 
to be considered. First, there is an 0( 1) outer wavgmmber region, 
where the stabilizing effects of surface tension, viscosity, and gravity 
are all relatively weak. Second, there is a far outer scale k ~k f /t , 
where surface tension and/or viscosity are strongly stabilizing and 
gravitational effects are, to a first approximation, negligible. Finally, 
we have an inner scale or*** 2 , where gravity is the dominant 

stabilizing effect (the first scale is valid for normal gravity, the latter 
for the reduced gravity regime) and where viscosity and surface- 
tension effects are absent at leading order. In each of these regimes, 
the cellular stability boundary, obtained by seeking solutions of 
the dispersion relation for which ia> is identically zero, is given, 
respectively, by 




i 


a;w ■ 


p_ 
' 2 


■ 0 ( 1 ) 


[P*(P*« -*»)/“»• Fr 
Ap -*)/*. Fr'~Ob) (35) 

WPli+WrWPWiLj fr 

” p + 4p*p(l +p'Pkf) - (1 - RKp'y +2p.'Pr~' 

7?=(l+4p* 2 rt}) 1 (36) 

where P a Fr/. Matching these solutions to one another then leads 
to the composite stability boundary 

a; ( c) - -p* 


2p* p' P [ l+€k(p*y +2p*P-+2p*P)] 


Although the preceding problem is linear in the coefficients 
hj-h io, which caiy'thus/be eliminated to give a single equa- 
tion for io), the resulting dispersion relation is quite long and 
highly nonlinear. Explicit results may be obtained for certain spe- 
cial cases, including the original problems considered by Landau 
(A p = Ae=Fr/=M = 0) and Levich 2 (A p - A© = ja = y = 0), as 
well as a particular case (A© = /x = Pn = 0) of die generalized 
model described earlier. 5 * 13 To obtain more general results, it is 
possible to exploit the smallness of certain parameters and to seek 
asymptotic solutions for the neutral stability boundaries. In partic- 
ular, realistic limits to exploit include the smallness of the gas-to- 
liquid density and viscosity ratios p and /a, and, in the microgravity 
regime, Ff~ , Pursuing this approach, tractable asymptotic results 
have so far been obtained for A© = 0 (Refs. 5 and 6) and for the 
inviscid problem when A© is nonzero (Ref. 7). The present work 
essentially completes the asymptotic analysis of the dispersion rela- 
tion embodied in Eqs. (27-34) By extending the last of these studies 
to the fully viscous case. 

Parameter Scalings and Asymptotic Analysis 
of the Dispersion Relation 

Focusing on the realistic regime p <SC 1 (typical values are on the 
order of 10' 3 or 10' 4 ), we formally introduce a bookkeeping param- 
eter € 1 and introduce the reasonable scalings p = p*£, p = p**, 

Pn~0( 1), and either Fr‘ l or where Fr- l ~0(O cor- 
responds to the case of greatly reduced gravity. In this parameter 
regime, the appropriate scaling for A p to describe the neutral stabil- 
ity region is A p = A* * (Refs. 5 and 6), whereas the appropriate scale 
that describes the main effects of thermal coupling turns out to be 
A© = A* € l/ * (Ref. 7). Based on this scaling, we note that the ratio 
Ae/A ~ ~ 3 / 4 ) 1, as might be expected based on an overall 

Arrhenius reaction-rate dependence on temperature. 1 

Based on our earlier analyses, the scalings introduced induce a set 
of corresponding regimes for the wave number k (and the complex 
frequency ito) in the dispersion relation determined by Eqs. (27- 


.4 p 

£ 

2k 


*P(l+«*p‘P)-(p*)' + 2p*/ , )[l-(l+4p ,J f >1 « 2 * J )*] ■ 
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«Y 


(37) 


)# 

A 


as shown in Fig. 1 . Clearly, the stable region lies below A* = -p* /2 
(negative values of A, over certain pressure ranges are characteristic 
of a number of HAN-based liquid propellants^ with the location 
of the minimum in the cellular boundary increasing to less negative 
values of A„ with increasing values jf the stabilizin g parameters 
yjP p r . ,i\ and * (or g‘) In Fig. I She curves arc drawn forthcv 
case e = 0.04. p* = 1 .0 j = 6.0^* S 2.0. Tlic solid curves corn- 1 
spond to (lie inviscid lifiiit (P = 0) wiUi nonzero surface tension 
(y = 2.5). The dash-dot curves correspond to nonzero surface ten- 
sion (y = 2.5) and liquid viscosity {P = l 0), but zero gas-phase vis- 
cosity (/a* P = 0). The dash-dot-dot curves differ from the dash-<lot 
curves by the addition of gas-phase viscosity (/a* P = 1*0) are 
similar to the dash-dot-dot-dot curves, where the latter correspond 
to larger viscosities {P = /t*P = 2.0). The dash-dot-dot-dot-dot 
curves correspond to a viscous case (P = p* P = 1 .0), but with zero 
surface tension ^fihe£ainp^ sets of curves cor- 

responding to tfw normal and reduced gravity cases, it is clear that 
the critical wave number for instability becomes small in the latter 
regime. That is, cellular hydrodynamic instability becomes a long- 
wave instability in the limit of small gravitational acceleration. Fur- 
ther discussion of this stability boundary, and its relationship to jthe 
original Landau i /Levich 2 predictions, is given in Ref. 5. 

Considering the pulsating stability boundary (in the limit 
0), which is obtained by seeking solutions of the dispersion 
relation for which only the real part of i d) vanishes, it is found Uuit 
the corresponding expressions in the inner and, num er 

regions are given by . 

a;~-p*. a;~-pXi;+2/>*)1 

respectively. In tliis case, it is clear dial the outer solution is, in fact, 
the comp osite solution, which lies below the cellular boundaries an d 

A) C(/ vu 2JP ^JA A 




tnishes, it is found 8 that! I 
and. outer wave number 

■■ ’ • i : \ r :,v.i 
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5 HYDRODYNAMIC STABILITY BOUNDARIES (p « 1) 

1 Viscous Case (P > 0) 
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Viscous Case (P > 0) 
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Fig. 1 Asymptotic representation of the cellular hydrodynamic neutral 
stability boundaries. 

recedes to negative values of A p that are larger in magnitude than 
C?(e) as k becomes large (Fig. 2). Clearly, this stability boundary 
is more sensitive to the stabilizing effects of the liquid viscosity 
parameter P than is the cellular boundary, having a leading-order 
stabilizing effect for 0(1) wav^jbmber disturbances in this case. In 
the limit P -► 0, the pulsating boundary collapses to the straight line 
A* = -p*, that is, A * = - 1 in Figs. 1 and 2, but even in that limit, 
there is a region of stability corresponding to values of A* p greater 
than -p* and less than the minimum in the cellular boundary, which 
is greater than -p*( 2. However, if one now considers the effects of 
a nonzero temperature sensitivity in the inviscid limit P~ 0, then, 
for A 0 ~ C?(€ 1/4 ), the pulsating boundary possesses a turning point 
such that the stability region disappears for sufficiently small wav^J) 
number perturbations. 7 This is shown in Fig. 3, which indicates that 
the pulsating boundary then frames the stable region except along 
the upper branch that asymptotes to the previous cellular boundary 
as k becomes large in the outer wa\^umber region. The evolution 
from a stability diagram that indicates a stable region delineated by 
distinct pulsating and cellular hydrodynamic stability boundaries 
to the pulsating-dominated one shown in Fig. 3 can be shown to 
occur in the parameter regime A & ~ C?(t 1/2 ), which, based on the 
estimate A e /A p ~ C?e“ 1/7 > 30 (i.e., of the same order as a typi- 
cal nondimensional activation energy), appears to be attainable for 
many types of liquid propellants. We now extend the analysis that 
produced the fully developed pulsating boundary shown in Fig. 3 to 
the viscous case in which both P and p* are allowed to be nonzero. 

Owing to the complexity of the fully viscous problem, we analyze 
Eqs. (27-34) directly by seeking appropriate expansions for the 
complex frequency i co and the coefficients b t . This differs from our 
approach in the inviscid limit where it was feasible to first eliminate 
; the bt to obtain a singleimplicit equation for i<o alone. We first 
consider the 0(1) wavAmiber region and, based on our earlier 
analyses, seek an expansion for the dispersion relation io>(k) in this 
region in the form • '* ** 

io> ** € + ♦ ■ • ) (39) 


Fig. 2 Asymptotic representation of the pulsating hydrodynamic sta- / j 

bility boundary for the viscous case (P > 0).., 7#e CC/ft <>• (+**#&**%' 

Introducing the already defined parameter scalings, die quantities |i 

p y q, /*, and s defined following Eq. (26) are expanded as % ij 

p^p,_$>«~* + Po + P ( })«* +■" 


9 ~ + <?o + -" 


where 




5 ~ s 


*!»(«*); 


(40) 


= 2p 0 = [l + wi/('"o) i ] 

8p,j> = (/<uo) _ l [l + 4 k l + 4ia)j - (iwi) 1 / /a*] 

= (iW^L 2P< t° = [‘ + 

'"(J) = P'- . r <J) = 

rj = -iu>ip* - (p* Pk) 1 ’i 

i 1 ’• ’ ' ‘ , 

Finally, the bi are conservatively postulated to have the expansions 
s , . t ;-".- 1 +»!"*’<-* ; i - 1 . 1 * 

i. : •; ; ; : * \>\ . - 

f-li i tn\ n . (-1) i f ^ n in , /A 1 \ ‘ : W* 


:\A r kV 
' \ 


our earlier analyses of more specialized cases, v! i« 

Substituting the preceding expansions into Eqs. (27—34) and; > <5 -.! 
equaUng coefficients of like powers of wc 0^®“ *he. leading- ^ iVJ 
order conditions ■ %: v • : *. ; It ■ 
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A% = .5, p’ =y= I, « = • 005 



ik^'+ik^+q^b^' 

= ,^r i, +‘^‘ i) +^i)"i' ,, = o < 42) 

b<-‘> = -Jf/p*, i'r" = -(i«X>)V* 

fc H>= t ;-i > =o ( 43 > 

+ fc«-’> - 2fc'- 1) = 0, b ( i n = -C*/p’)(l + a;/ P *) 

(44) 

where the last of these was obtained from the leading-order differ- 
ence of the first and second of Eqs. (29) using the last of Eqs. (43). 
From Eqs. (43) and (44) we/thus/ conclude that the leading-order 
dispersion relation is given by 

(icoo) 2 -(A/P*) 2 (2AJ+P*) < 45 > 

which is the same result as that obtained in the inviscid case 
when Ac =0. In particular, Eq. (45) implies that (<a>o) ^0 for 
A* ^ _p* /2 y which recovers the leading-order cellular boundary 
(3^) for 0(1) wavenumbers, but gives no definitive information re- 
garding stability for A* < — p*/ 2 because ico o is purely imaginary 
in that region. That is, the stability of the basic solution in the latter 
region is determined by the real parts of higher-order coefficients in 
the expansion (39) for ia> , although lm[i a*} ^ 0 implies that distur- 


conclude that b { { m) =ia> od - A' P /p'), However, jK««i =0 im- 
plies the need to continue the analysis at the next order to determine 
ia> 2 . Proceeding in this fashion, we obtain from the earlier results and 
Eqs. (29-34) atthis next higher order anew set of conditions given by 

• - (k/iwo)b\~^ —2i<Pi = i(l - A*/ p% ] bf = 0 
’ '' (49) 

(50) 




X.C0) . 


- 2b { y y = - 1 + A'jp*) 

b^ -{K/P')b^ = (si) 

-2i^ = 2 p k \ \ ^ - - <J) 


H 

■’i 

i 


stability boundary. 

At the next order in the analysis of Eqs. (27-34), 
second set of conditions given by 

we obtain a 


+ ik 2 b { f^ * >=0 

(46) 

ik^ + ik^+r^+r, h bt l) =0 
iwi = b x = b 2 = = °s u i 

(47) 

= b? = C = = 0 

(48) 


= b\ * 

(52) 

where Eq. (5 1) was actually obtained from the next higher-order dif- 
ference of Eqs. (29) and the second of Eqs. (49) was obtained from 

die sum of Etyf. (31) multiplied by it, and E</. (32) multiplied bjr J / 
Equations (49-52) constitute a closed system for i>, ,u % 

b (l/4> and ia>i. Eliminating tile first three of these in favor of die last 
and using the result (45) for io*. we finally obtain the dispersion 
1 relation for icuj as 
iaj, = —IPk 2 + t(4*/ p' — l) 

x [a*/p* + 1 -t- P* - **-l/le(2A* / p* -h l) *] (53) 

Stability in die region A* < -/>*/ 2 below die cellular boundary 
is determined by the real part of ic^j. In that region, the principal 
value of the complex factor in Eq. (53) may be written as ^ApJ 
♦ + l)-3/4--[_(AVp*-H)]“ 3/ V 3<,r/4 , and thus the neutral 
stability condition 'Re{i<uj} = 0 leads to an implicit equadon 
for the (pulsating) neutral stability boundary A* (k; A e , P). In 
terms of the new pressure sensitivity parameter b defined by 
a* — _(«* /2)(1 + h), where b represents die negative deviation, in 
units of pV2/romlhe cellular boundary A* p = -p*A this boundary 
is given by i 4 > * ^ ’ 


1 


b* (3 + b)“ J l(3 + b)(l — i) + 8 Pkf =«*/ 


where the last of Eqs. (48) was deduced from the next-order dif- 
ference of Eqs. (29). Finally, from the sum of the first of Eqs. (28) 
muldplied by iki and the second of Eqs. (28) multiplied by lij, we 


• (54) 

L • 

In the limit k -+ oo. it is clear dial tha^re two solutions of Eq (54) 
given by b = 0, dial is. A* = -p 'll, and b ~ - 1 + 2( 1 + 2 Pk) , 
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that is, A*/p * (1 + 2Pk) l/2 . Thus, the pulsating boundary is 

clearly multivalued, as in the inviscid case (Fig. 3), with one branch 
approaching the cellular boundary and the other branch approaching 
the pulsating boundary for Ae — 0 (Fig. 2) in the limit of large k. 
More generally, Eq. (54) may be rewritten as a cubic equation for 
the inverse relation k(b) as 


64P 2 * 3 + 16(3 + b)( 1 - b)Pk 2 + (3 + b ) 2 { i 
-ai(3 + £) 2 £>~* =0 


-b) 2 k 


(55) 




which is clearly seen to collapse to the previous inviscid result 7 in 
the limit P-+0. For arbitrary P t typical plots of k(b) are shown 
in Figs. 4a~4d, which, when rotated -90 deg so that the k axis is 
horizontal, is readily interpreted in the context^of Figs. 1—3, where 
the lines A* - -p */ 2 and -p* correspond to b = 0 and A J, respec- 
tively. It is clear that thesexurves asymptote to the lines b = 0 and 


.V uveiy. U la vac<u umi ■ — — j — , 

r q. 2(1 + 2 Pk ) 1/2 as k oo, where the latter corresponds to the 

^viscous pulsating boundary in the limit A * e — ► 0. They cross the 
* line b - 1, which corresponds to the inviscid pulsating boundary in ■, 
jfrcrf .to cpeeified'limit, at k* = a 1/2 /4F 2 . That the pulsating boundary 
A becomes C shaped (in the rotated frame of reference) for A e > 0 
implies that steady, planar burning is intrinsically unstable for suffi- 
ciently small wave numbers. In addition, because the portion within 


the C-shapcd curve is lire stable region, any crossing of the C-shaped 
boundary from the stable to the unstable region corresponds to the 
onset of a pulsating instability. As A© increases, the turning point of 
the C-shaped pulsating boundary shifts to larger values of k. On the 
other hand, as A* e becomes small, the turning point shifts to smal^ 
values of k such that this point eventually leaves the (2(1) wav^/ 
number region for which Eq. (54) is valid. Indeed, it turns out that 
the transition to separated pulsating and cellular branches occurs as 
A e decreases through <2<€ ,/2 ) values for intermediate (2(e wave 
numbers. 7 Thus, as A% becomes small, the original pulsating and 
cellular boundaries are recovered in the (2(1) wave number regime, 
but as AJj becomes large, the original cellular boundary lies within 
the unstable region for (2(1) wave numbers, and the basic solution 
becomes intrinsically unstable to oscillatory disturbances. 

Composite Neutral Stability Boundary 
A composite asymptotic solution for the neutral stability bound- 
ary in the regime A 0 - <2(€ 1/4 ) is thus obtained by matching die cel- 
lular and pulsating boundaries in the far outer wave number regime, 
where the former is given by Eq. (36) and the latter by the second 
of Eqs. (38).widi die appropriate solution branch of Eq. (54) in the 
(2l>^wav^imber region. In particular, reverting back to the pa- 
ram^tfer A* ?we denote the two solution branches of Eq. (54), which 
correspond to die portions of Fig. 4 dial lie to the left and to the 
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Fig. 5 Composite pulsating/cellular hydrodynamic stability boundary for Aq ~ 0(e 1 ^ 4 ) in the general viscous case* ? 


right of the turning-point minimum, by A* p ^°'*\k) and A* p +*(£)• 
\*where the superscript o jtfTas before, the outer, or 0(1), wavtomm- ^ 
* ber region and the superscripts u and / ^the upper and rower 
(rotate Fig. 4 by —90 deg) portions of the double-valued pulsating 
boundary A*(k ). Along the upper branch, A p {Q,u) -p*/2 t that 
is £ — ► 0, as k — ► oo, which can be matched with Eq. (36) because 
A*V> - P 72 as k f -* 0. Similarly, A;* 0 -+ -p*(l + 2Pk ) l/2 
[i.e., b-+ -1 + 2(1 + 2 Pk)' /2 ] as *-► oo, which clearly matches 
the viscous pulsating boundary given by the second of Eqs. (38) 
in the far outer wave number region. As a result, a leading-order 
composite stability boundary spanning both the outer and far outer 
wav^umber regions is given^y 


li-y 

\ v * 




with the combustion front, realistically possesses both a pressure and 
temperature sensitivity. In the present work, the fully viscous case 
was considered, thereby generalizing previous analyses in which 
either the viscosity of the fluid and/or the temperature sensitivity 
of the reaction rale was neglected. As in these preceding studies, 
the smallness of the gas-lo-liquid density ratio was used to define 
a small parameter that allowed an asymptotic treatment of a rather 
complex dispersion relation. Specifically, it was again shown that 
in addition to the classical Landau, or cellular, stability boundary, 
there exists a pulsating hydrodynamic stability boundary as well. For 
sufficiently small values of the temperature-sensidvity parameter, 
there is a stable region between these two boundaries corresponding 


p* 2p'p'P[l -M*(pV +2/x»P +2p»F)] 

V ' *(*> - y + 4 M .p(j + P 'Pek) - (l - [1 )(/»*>' + 2 m* P) 

A'M{k) 


for AySU* c , where A p c iyihe turning point calculated from 
Eq. (MjanS the second tenfi in the top expression has been ex- 
pressed in terms of the out^av^umber variable k. 

The composite stability boundaiyls shown in Fig. 5. Based on the 
preceding construction, the lower branch of Eq. (56) is a pulsating 
boundary for all wave numbers, whereas the upper branch transitions 
from a pulsating boundary for 0(1) wave numbers to a cellular 
boundary for 0{e~ l ) wave numbers. Indeed, from Eq. (45), the size 
of the upper region of oscillatory instability, which is bounded below 
by the upper branch of the pulsating stability boundary and above 
by the region of nonoscillatory instability beyond the outer cellular 
boundary A' ~-p*/2 for A' e =0, shrinks to zero as k becomes 
large on the 0(1) wav^gjumber scale. In this regime, the lack of a 
stable region for sufficiently small wave numbers/thus / /fmplie$ an 
intrinsic instability to long-wave pulsating perturbations. 

Conclusion 

. ' The present work further extends our recent formal treatment 
of hydrodynamic instability in liquid-propellant combustion. The 
analysis is based on a generalized Landau/Levich model in which 
the dynamic motion of the liquid/gas interface, assumed to coincide 


to a range of negative pressure sensitivities for which steady, planar 
burning is stable. 

As the pressure sensitivity decreases in magnitude, the cellular 
stability threshold is crossed, leading to classical Landau instabil- 
ity. Surface tension, viscosity (both liquid and gas), and gravity 
are all stabilizing effects with respect to this type of instability. 
However, only gravity stabilizes small wavg)umber disturbances, 
and thus Landau instability becomes a long- wave instability in the 
reduced-gravity limit. Alternatively, as the pressure-sensitivity pa- 
rameter increases in magnitude, the pulsating boundary is crossed, 
and liquid-prppellant combustion becomes unstable to oscillatory 
perturbations. This type of hydrodynamic instability is more sensi- 
tive to the stabilizing effects of (liquid) viscosity than is the cellular 
boundary, but the stabilizing influence of viscosity does not ex- 
tend to small wav^umber disturbances, and gravity turns out not ^ 
to have a significant effect on this type of hydrodynamic instabil- 
ity. Consequently, for sufficiently large values, of- the temperature-.^ 
sensitivity parameter, the pulsating boundary 'develops a turning J 
point and becomes C shaped. In this parameter regime, correspond- .* 
ing to ratios of the temperalure-to-pressure sensitivities of the order 
of 200-1000, steady, planar combustion is intrinsically unstable to 
nonslcady long-wave perturbations. In that casc v the pulsating form 
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of hydrodynamic instability is predicted to dominate, leading to 
large unsteady cells along the burning liquid/gas interface. 
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